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Abstract — Cadambe and Jafar (CJ) alignment strategy for the 
if -user scalar frequency-selective fading Gaussian channel, with 
encoding over blocks of 2n + 1 random channel coefficients (sub- 
carriers) is considered. The linear zero-forcing (LZF) strategy 
is compared with a novel approach based on lattice alignment 
and lattice decoding (LD). Despite both LZF and LD achieve 
the same degrees of freedom, it is shown that LD can achieve 
very significant improvements in terms of error rates at practical 
SNRs with respect to the conventional LZF proposed in the 
literature. We also show that these gains are realized provided 
that channel gains are controlled to be near constant, for example, 
by means of power control and opportunistic carrier and user 
selection strategies. In presence of relatively-small variations in 
the normalized channel coefficient amplitudes, CJ alignment 
strategy yields very disappointing results at finite SNRs, and 
the gain of LD over ZLF significantly reduces. In light of 
these results, the practical applicability of CJ alignment scheme 
remains questionable, in particular for Rayleigh fading channels, 
where channel inversion power control yields to unbounded 
average transmit power. 

I. Introduction 

"Everyone gets half the cake" is the surprising promise of 
interference alignment, as introduced by Cadambe and Jafar 
[JTJ for a A'-user Gaussian interference channel with ran- 
dom coefficients. Interference alignment is a linear precoding 
strategy that forces the interfering signals at each receiver 
k to span a subspace Ik of the receiver signal space, such 
that the desired signal can be transmitted in a subspace Sf. 
with If. U Sk = {0}. Therefore, each receiver k can remove 
interference completely by a linear zero-forcing projection on 
the orthogonal complement I^r. With enough richness in the 
channel coefficients, in the form of jointly distributed random 
variables drawn from a continuous distribution (e.g., as arising 
from time and/or frequency selective fading channels) (TJ 
shows that by encoding over a large block of N channel 
uses, the limit of K/2 degrees of freedom is achievable, 
i.e., lini7v->oo dim(Sk)/N = 1/2 for all k. With such a 
strategy in place, the entire group of interfering transmitters 
at each destination appears collectively as a single source of 
interference, consuming only half the total degrees of freedom 
(dimensions) and leaving another half available at each user. 

In this paper, we put this idea to test by investigating practi- 
cal decoding strategies for interference alignment in practical 
signal-to-noise ratios (SNR). We focus on the parallel channel 
single-antenna scenario introduced in JT1 for i.i.d. random 
subchannel coefficients, e.g., OFDM with frequency-selective 
fading. We find that despite the promise of the degrees-of- 



freedom analysis, the quality of the effective channel with 
interference aiglnment and linear zero-forcing (LZF) interfer- 
ence removal is generally not quite acceptable at practical 
finite SNRs. This is especially true for channel coefficients 
with non-constant amplitudes (as in a frequency-selective 
channel). 

The major limiting factor in finite-SNR performance of 
Cadambe and Jafar (CJ) interference alignment pertains to 
channel inversion operations at the transmit side precoding, 
and zero-forcing interference cancelation at the receive side. 
In order to improve upon the CJ strategy, in this work we 
consider discrete alignment of signal sets in addition to the 
CJ alignment of signal subspaces. The core idea is to obtain 
an equivalent MIMO channel suitable to more efficient non- 
linear decoding strategies. In order to do so, it is essential that 
the superposition of interference signals not only span a low- 
dimensional subspace, but also appear as points in a discrete 
lattice, closed with respect to addition. The advantage of 
discrete lattice alignment is that it allows to recast the channel 
observed by each receiver as a MIMO channel where the 
desired signals occupy approximately half of the signal space 
dimensions, and the sum of interfering signals spans the other 
half. This allows to apply well-understood Lattice Decoding 
(LD) strategies such as sphere decoding and variations thereof 
(e.g., (4)-|(6l), in order to decode the intended signal along 
with the superposition of interfering signals. 

We show that by exploiting the discrete nature of the 
transmitted signals, this nonlinear decoding approach yields 
a surprising improvement compared to linear alignment with 
zero forcing, provided that the dynamic variations of channel 
coefficient amplitudes are controlled. Channel dynamics have 
degrading effects on the performance largely due to the 
transmit-side inversive precoding. However, suitable carrier 
pairing in time and frequency dimensions, and opportunistic 
user selection strategies like those in JT) could be used as 
powerful means to achieve near constant channel amplitudes. 

The rest of the paper is organized as follows. Section [II] 
gives a summary of the CJ interference alignment scheme for 
single-antenna, parallel fading channels. The idea of discrete 
alignment and the MIMO interpretation of interference chan- 



nel with alignment is introduced Section III and Section |TVj 
presents performance results and simulations. Section |V] con- 
cludes the paper with a few final remarks. 



II. Linear Subspace Alignment 

We focus on the K = 3 user case, with precoding block 
length N — 2n + 1 for some integer n > 1, as in pi. A 
channel use of this channel is given by 



Y* = Hl x X\ + H\ 2 X\ 



^2 — H 2i X 1 + H 22 X 
Y 3 = H 3l X 1 + H 32 X. 
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where {Hjj} denote the channel coefficient, and {Y*, JT|, Z\\ 
represent the received symbols, the transmitted symbols, and 
the noise samples at channel use t, respectively, for i,j = 
1, 2, 3. The CJ strategy is a precoding scheme across a block 
of dimensions (e.g., subcarriers in an OFDM channel). In 
a block of length 2n + 1 channel uses, user 1 encodes a 
vector Xi of length n + 1 symbols using a precoding matrix 
V 1 eC (2n+1)x( " +1) , while users 2 and 3 encode X? and X3 
of length n using precoding matrices V 2 , V 3 6 02n+i)xn^ 
respectively. Notice that the encoding scheme is not com- 
pletely symmetric. However, symmetry can be achieved on 
average, by a rotating scheduling among users over successive 
blocks. 

The equivalent block channel can be described as 



Yi — Hn V1X1 + H12V2X2 + H13V3X3 
Y2 = H21V1X1 + H22V2X2 + H23V3X3 
Y 3 = H31 V1X1 + H32V2X2 + H33V3X3 



Zi (la) 
Z 2 (lb) 
Z 3) (lc) 



where Yj, Zj represent the received and noise vectors of length 
2n + 1, Xj represent the users' data vectors, and denote 
the diagonal channel matrix of size (2n + 1) x (2n + 1), for 
i,j = 1,2,3. 

The precoding matrices Vi,V2,V3 are designed such 
that the interfering signals occupy a common subspace. This 
alignment is achieved in [1] through the following design: 

(2n+l)xl 
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H 32 i H 31 V 1 P 2 
H 23 1 H 2 iV 1 P 3 
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It is fairly straightforward to check that with the above 
precoding design, we have 



H12V2 = H13V3, 

H23V3 = H^iVxP^, 

H32V2 = H 3 iViP 2 , 



(4) 



where (a) follows since by the structure of T and the shift 
property of P2 and P3, we have V1P2 = TV1P3. As a 



consequence, the two interfering signals received at user i are 
confined to a common subspace spanned by the columns of 
H12V2, H21 Vi, and H31V1 for i = 1,2, 3, respectively. This 
is readily seen if we rewrite (??) as 

Yi = HnViXi + Hi 2 V 2 (X 2 + X 3 ) + Zj (5a) 
Y 2 = H 2 2V 2 X 2 + H 2 iVi(X a + P3X3) + Z 2 (5b) 
Y 3 = H33V3X3 + H 31 Vi(X a + P 2 X 2 ) + Z 3 . (5c) 

The aligned structure of the interference vectors allows the 
decoders to remove interference by LZF, i.e., by projecting 
the received vector onto the orthogonal complement of the 
interference subspace. 

Remark 1. For successful zero forcing, we also need the 
interference subspace to be linearly independent of the signal 
subspace, i.e., [ HuVi Hi 2 V 2 ], [ H 22 V 2 H 2 iVi ], 
and [ H33V3 H31V1 ] in |5]) are full rank. In [1], it is 
shown that this rank constraint is almost surely satisfied for 
channel coefficients drawn from a non-degenerate continuous 
distribution (i.e., a distribution for which no coefficient is a de- 
terministic function of other coefficients). Channel coefficients 
are further constrained in nlyto satisfy a < \Hfj\ < b, for 
some < a < b < 00. (see|7] Section II].) This magnitude 
constraint, in particular, rules out direct application of this 
alignment strategy in a Rayleigh fading environment, for which 
channel inversion yields unbounded average transmit power. 
However, this strictly-bounded constraint can be enforced, in 
practice, by power control and opportunistic user selection. 
For example, we may assume that each subcarrier t in an 
OFDM system is pre-multiplied by a power control function 
\fPl such that for all t and all users i, \H t ji\fV\\ is bounded, 
for all j. Notice also that simple channel inversion does 
not accomplishes constant channel amplitudes, since a single 
transmitter must equalize the power of all its K outgoing 
links simultaneously. However, we may imagine some form 
of opportunistic subcarrier pairing across multiple users ^7^, 
such that groups of 2n + 1 subcarriers are chosen in order 
to have roughly the same amplitudes. By doing so, we have 
that j if*! I ~ \H\ 2 1 ~ l-^isl f or a M i — 1)2,3, and the power 
control command V\ is used to equalize amplitudes across the 
three transmitters. 



III. Discrete Alignment with Nonlinear Decoding 

Qualitatively, the key element in interference alignment 
is constricting the "expansion" resulting from linear super- 
position by confining the interfering signals to a common 
subspace. For discrete signals, however, we can prevent this 
"expansion" if the codewords form a discrete additive group. 
For example, consider a lattice A = {Gx|x € Z n }, where G 
is a full rank generator matrix. Summation of two codewords 
Xi , X2 G A gives another codeword xi + X2 G A. 

Consider for example the following interference channel 



Y l = g 1 X l + (X 2 + X z ) + Z± 
Y 2 = g 2 X 2 + (X 1 + X 3 ) + Z 2 
Y 3 = g 3 X 3 + (X 1 +X 2 ) + Z 3 , 



(6) 



where X\ , X 2 , X 3 are chosen from a rectangular QAM con- 
stellation C, given as a subset of the complex integer lattice 
Z[j]. Receiver i attempts to recover the desired symbol Xi 
by treating Xj + Xk (for j, k ^ i) as points in an expanded 
constellation, obtained as the sum set C = {x G C : x = y + 
z, (y, z) G C 2 }. In general, if symbols are uniform over C, the 
resulting extended constellation C is used with non-uniform 
probability; however, for the time being, we neglect this fact 
and consider maximum-likelihood (ML) decoding (assuming 
uniform prior probability for Xi G C and Xj + X^ G C). A 
more general one-dimensional strategy is proposed in [8| and 
is shown to achieve the maximum asymptotic degrees of free- 
dom for a scalar interference channel with deterministic (fixed) 



coefficients. See also [9], [10] where alignment strategies using 
lattice codes are presented for an interference channel with 
(up-to-a-scaling-factor) rational channel coefficients. 

A. Lattice Decoding 

We can combine the CJ "linear space" alignment strategy 
with the above (algebraic) lattice alignment idea. Let the 
data symbol vectors Xi,X2,X3 take on values in Cartesian 
product subsets of the cubic lattice with uniform probability, 
i.e., we let Xj G C n+ \ X 2 G C n and X 3 G C n , with C C Z\j]. 
Therefore, the interference terms X 2 + X 3 , Xi + P3X3 and 
Xi + P2X2 in |5]l are vectors drawn from the complex cubic 
lattices Z"[j] and Z n+1 [j], respectively. 



Rewriting |5]l as 
Yi = [ H 11 V 1 H 12 V 2 ] 

:= G1X1 + Zi 
Y 2 = [ H 22 V 2 H 2 iVi ] 

;— G 2 X 2 Z 2 
Y 3 = [ H 33 V 3 H31V1 ] 

:= GaXa + Z 3 , 
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user 1 can decode Xi along with interference sum X 2 + X 3 . 
Similarly, user 2 and user 3 decode their intended messages X 2 
and X 3 along with a shifted sum of interference Xi + P 3 X 3 
and Xi + P 2 X 2 , respectively. 

For Gaussian noise, ML decoding of Xi, X 2 , X 3 amounts 
to minimum square Euclidean distance decoding. For example, 
for user 1 this becomes 



xi = arg mm 

xee«+ 1 x(C)" 



lYi-Gixll 



Then, the decoder retrieve the first n + 1 components of Xi 
as the decision on the desired symbols Xi. This minimum 
distance decoding is identical to the MIMO decoding in a 
(2n + 1) x (2n + 1) AWGN MIMO channel, for which a 
vast amount of research exists. In particular, since Gi, G 2 , G 3 
are full rank (at least with random channel coefficients), an 
efficient strategy for minimum distance search consists of 
ignoring the finite constellation constraints and search over the 



entire cubic lattice x G Z ™ + [7], using well-known closest 
lattice point search algorithms, commonly known as "sphere 
decoding" (see (4j-||6| and references therein). 

B. Linear vs Nonlinear Interference Cancellation 
Consider user 1: 



Yi — HnViXi + Hi 2 V 2 (X 2 + X 3 ) 
= GnXi + Gi 2 (X 2 + X 3 ) + Zi. 



■Zi 



(8) 



Consider ML decoding by treating the interference X 2 + X 3 
as arbitrary unknown vectors. This is equivalent to assume 
Gi 2 (X 2 + X 3 ) = S, where S is an arbitrary unknown vector 
with the constraint S € Span(Gi 2 ), i.e., S is any vector in 
the linear space spanned by the columns of Gi 2 . 

Under such assumption, the standard approach consists of 
the so-called Generalized Likelihood Ratio Test (GLRT): 

xi = arg max max p Yl \x u s(Yi\x.i,s), (9) 

xieC"+! se Span(G 12 ) 

which is equivalent to: 

xi = arg min min || Yi — GnXi — s|| 2 . (10) 

x 1 ec"+ 1 se Span(G 12 ) 

Writing s = Gi 2 u, with u G C", the inner minimization is a 
standard Least-Squares problem, whose solution is given by 



u = fGt 2 G 12 ) G f 12 (Yi - Gi lXl ) 



(11) 



Replacing this into the objective function of the outer mini- 
mization, we obtain 



xi 



arg min \\Pf 2 (Yj. - G a xi) , (12) 



where 



V 1 - 



I 



G12 ( G} 2 Gi 2 



'12 



is the orthogonal projector onto the orthogonal complement of 
the interference space Span(Gi 2 ). This projection corresponds 
to the familiar linear ZF receiver assumed in Cadambe-Jafar. 
Therefore, not surprisingly, the linear ZF receiver is the GLRT 
receiver assuming an arbitrary interference vector whose only 
constraint is to lie in a certain subspace. 

Now, we assume that the symbol vectors Xi,X 2 ,X 3 take 
on values in some discrete sets with uniform probability. The 
MAP decoder for X x takes on the form 



xi = arg max 



arg max 



(X! =x 1 |Y 1 ) 
Pyx\x x ( y iI x i) 



*iee«+i £xiec™+i p Yl \ Xl (Yi |xi) 



arg max 

x 1 6C»+ 1 



Sx, x, exp (- ||Yi - Gnxi - G i2 (x 2 +x 3 )|| 



E x ' ieC " + i^i|^i( Y l| x i) 



Taking the log and neglecting irrelevant terms, we have 



xi 



are; max 
x 1 ec"+ 1 



log 



E exp 

\X2,X 3 



lYi-GiiXi 




At this point, we notice that for sufficiently large SNR the log- 
sum of exponential terms is dominated by the largest term, 
which corresponds to the smallest distance. Therefore, we 
arrive at the approximated MAP decoder 

2 



Xl 



arg mm mm 



||Yi-Ghx 1 -Gi 2 (x2 + x 3 ) 



which, again not surprisingly, coincides with the advocated 
LD approach when X2 + X3 is a lattice point. 

We conclude that the LZF and the LD approaches corre- 
spond to ML decoding under arbitrary interference constrained 
into some linear subspace, and approximated MAP decoding 
taking into account the true discrete nature of the interference. 

IV. Performance 

It is true that the full rank criteria is almost-surely satisfied 
for channel coefficients drawn randomly from a continuous 
distribution. However, in addition to rank criterion, the overall 
performance of IA with LZF highly depends on the orthogo- 
nality of the channel coefficient matrices Gi, G2, G3. Notice 
that the channel matrices are ill-conditioned if two elements of 
the diagonal matrix T are close in value (and the rank criteria 
breaks if two elements of T are exactly the same). Further, the 
channel matrices are also ill-conditioned if the absolute value 
of any of the channel coefficients deviates from one, especially 
for larger n. The reason becomes clear by inspecting the 
exponential structure of the precoding matrices in |2} and the 
diagonal channel ratio matrix T in ([3]); if any of the diagonal 
elements of T have an absolute value lower (or higher) than 1, 
the corresponding row in Vi exponentially tends to zero (or 
infinity) for large n. Thus, we expect a poor performance from 
linear alignment with LZF for dynamic channel amplitudes. 

Consider first a random channel where channel coefficients 
all have magnitude 1 with a random uniform phase, i.e., Hh = 
exp{j>*J for i.i.d. 0*, ~W(0,27r), k,l = 1,2,3, 3 = sf-L. 
For this setup, Fig. [T] shows a simulation comparison between 
filtering out the interference subspace (LZF), and a nonlinear 
strategy using Schnorr-Euchner sphere decoding strategy (LD) 
of J5J. In this figure, the vertical access represents symbol 
error rate, and the horizontal axis represents average SNR. 
Here, SNR is defined as the ratio of average transmit power 
to noise, and is computed by averaging the realized transmit 
power (including the channel inversion) over a set of con- 
secutive blocks for fixed constellation energy and received 
noise variance. We compare the linear alignment strategy 
with LZF, and the discrete alignment scheme with LD, for 
blocklengths N = 2n + 1 = 21 and N = 2n + 1 = 11. In 
all simulations, the source symbols are drawn from a 4-QAM 
(QPSK) constellation. 

Fig. [T] shows that LD obtains a significant improvement 
over LZF. As shown in this figure, even with fixed channel 
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Fig. 1 . Comparison between LZF and sphere lattice decoding (LD) for fixed 
channel amplitudes. The source constellation is 4-QAM. 
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Fig. 2. Performance severely degrades with dynamic truncated channel 
amplitudes. Here, channel amplitudes are truncated within [0.8,1.2] and 
[0.8, 1.4], The source constellation is 4-QAM. Blocklength is N = 11, i.e., 
n = 5. 



amplitudes, performance of LZF is not quite acceptable. For 
SNR values as high as 80 dB, the symbol error rate is about 
10 -3 for n = 5. However, discrete alignment with sphere 
decoding achieves a reasonable symbol error rate in the order 
of lO- 3 for SNRs around 20-25dB. 

Theoretically, larger blocklengths achieve higher degrees 
of freedom. However, for the purpose of decoding at finite 
SNRs, larger blocklengths result in a higher symbol error 
rate, because of the skewed structure of the equivalent MIMO 
channel. This can be seen by comparing the corresponding 
curves for both LZF and LD at n = 5 and n = 10 in Fig. [T] 
Recall that the structure of the equivalent MIMO channel 
matrix is exponential in n. 



Fig. [2] shows a similar comparison between linear alignment 
with LZF and discrete alignment with lattice decoding (LD), 
for truncated Gaussian coefficients with amplitudes within 
[0.8, 1.2] and [0.8, 1.4]. In this simulation, channel coefficients 
are drawn from a complex circularly-symmetric normal distri- 
bution, and subchannels with a magnitude outside the desired 
interval are discarded. In comparison to symbol-error rate 
performance with fixed channel amplitudes in Fig. [T] dynamic 
channel amplitudes significantly degrade the performance as 
shown in Fig. [2] In fact for larger block lengths, the advantage 
of LD over LZF quickly fades out with increasing channel 
amplitude dynamics. 

It is intuitively predictable that smaller channel amplitudes 
directly result in worse performance. However, Fig. [2] shows 
that larger channel amplitudes are also not suitable for in- 
terference alignment and result in a significant performance 
loss. As shown in Fig. [2] LZF performs very poorly even 
at SNRs as high as 35dB for dynamic channel amplitudes 
within [0.8, 1.4]. However, with discrete alignment combined 
with LD, we obtain a reasonable symbol error rate of 10 -2 
(which combined with a powerful outer code could yield an 
acceptable overall bit-error rate) at SNR 35 dB. A noticeable 
improvement is further observed when channel amplitudes 
are truncated within [0.8,1.2], which again, demonstrates the 
effect of channel dynamics on the overall performance. We 
conclude that the precoding strategy in |2]i and |3]l is highly 
susceptible to power control and channel amplitude dynamics, 
and a tight power control strategy and carrier pairing is 
required to achieve a reasonable performance. 

V. Conclusion 

While the asymptotic degrees-of-freedom analysis is a 
powerful tool in understanding sophisticated communication 
systems, it provides little intuition about performance at finite 
SNRs. In fact in many cases of interest, a one-dimensional 
(scalar) analysis suffices to compute asymptotic degrees of 
freedom, including the if -user interference channel [8], relay 
networks ( TT) , the two-user interference channel [12], and 
many-to-one and one-to-many interference channels (T3). 

In this paper, we examined the finite-SNR performance 
of linear subspace interference alignment strategy of [ 1 1 for 
frequency-selective interference channels with single transmit 
and receive antennas and i.i.d. continuous channel coefficients. 
A discrete alignment strategy is introduced that uses algebraic 
properties of lattices to transform the single-antenna inter- 
ference channel to an equivalent MIMO system, for which 
efficient decoding strategies are well understood. While linear 
subspace alignment with zero forcing shows a poor perfor- 
mance at finite SNRs, discrete alignment along with sphere 
decoding strategy significantly boosts the performance. 

We observed a particular performance sensitivity to channel 
amplitude dynamics. This suggests that to benefit from inter- 
ference alignment, we need to employ a tight power control 
mechanism, combined with user selection, and carrier pairing 
in time and frequency dimensions to obtain near-constant 
channel amplitudes. 



The ergodic interference alignment strategy of G) avoids 
exponential beamforming vectors and channel inverse opera- 
tions, thus, could be more suitable for larger channel amplitude 
dynamics, e.g., Rayleigh fading. In ergodic interference align- 
ment, the transmitting nodes, in a way, employ an extreme 
form of subcarrier pairing over the time horizon, essentially 
"waiting for the right channel" to align. As a result of this 
opportunistic approach, however, ergodic alignment suffers 
from an average decoding delay that grows exponentially as 
0(exp(BK 2 )), where B is the number of quantization bits 
used to represent the channel coefficients (which also grows 
with the operating SNR). On the other hand, the CJ alignment 
strategy ignores the inherent randomness of channel coeffi- 
cients and employs a deterministic precoding strategy that 
works for a fixed deterministic blocklength, though the penalty 
in power could be large. The next step to further improve 
the performance could be using opportunistic strategies to 
control channel variations to obtain new interference alignment 
strategies in between the CJ scheme and ergodic alignment. 
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